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Abstract 

We develop the quantization of unimodular gravity in the Plebanski and Ashtekar 
formulations and show that the quantum effective action defined by a formal path 
integral is unimodular. This means that the effective quantum geometry does not 
couple to terms in the expectation value of energy proportional to the metric tensor. 
The path integral takes the same form as is used to define spin foam models, with 
the additional constraint that the determinant of the four metric is constrained to be 
a constant by a gauge fixing term. This extends the results of [10 j to the Hilbert space 
and path integral of loop quantum gravity 

We review the proposal of Unruh, Wald and Sorkin- that the hamiltonian quan- 
tization yields quantum evolution in a physical time variable equal to elapsed four 
volume-and discuss how this may be carried out in loop quantum gravity. 
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1 Introduction 

The unimodular formulation of general relativity was first proposed by Einstein in 1919 
as an approach to the unification of gravity and matter|[T]] . It was studied by a number of 
authors in the 1980s and early 90s because of indications that it resolves two key problems 
in quantum gravity[2J-[9J. These are the cosmological constant problem and the problem 
of defining a physically meaningful time with which to measure evolution of quantum 
states in quantum cosmology, in the absence of a spatial boundary. This is the second of 
two papers which report results which support and clarify the sense in which unimodular 
quantum gravity solves these two problems. 

In the first of these papers [TTJ|, I constructed the constrained phase space quantization 
of a formulation of unimodular gravity due to Henneaux and Teitelboim[8J. I showed 
that the quantum effective action is a functional of the unimodular spacetime metric 
with determinant fixed to 

\Jdet{gy, v ) = e (1) 

where eo is a fixed nondynamical volume element. This means that the quantum effective 
equations of motion, which arise from varying the metric with ([1]) fixed have a symmetry 

T ab -> T' ab = T ab + g ab C (2) 

where C is a spacetime constant. This decouples the dynamics of the metric g from any 
contribution to the energy-momentum tensor, whether classical or quantum, of the form 
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of a constant times the spacetime metric. This means that the puzzle of why huge contri- 
butions to of this form coming from the zero point energy of the fields, expected to be 
of the form of Mp t are not sources of spacetime curvature. This is sometimes called the 
old cosmological constant problem. 

In this paper these results are rederived using the Ashtekar variables and the Ple- 
banski form of the action. Specifically, in sections 2 to 4 we develop the unimodular 
formulation of the Ashtekar variables first studied in HI. In section 5 we formally con- 
struct the path integral and show that the quantum effective equations of motion have 
the unimodular property ©. 

This sets up the possibility of explicitly realizing unimodular quantum gravity in the 
context of a spin foam model, which is a well defined and ultraviolet finite path integral 
quantization of general relativity. 

In section 6 we then seek to build on these results to address the problem of time 
in quantum cosmology and the related issues of defining physical observables and the 
physical inner product. Here we are following a suggestion made by Unruh[5|], Wald|6], 
Sorkin[7l and others, which can be summarized as follows. We define, following Hen- 
neaux and Teitleboim a three form a a ^ 7 which satisfies a field equation 

b = (da)* = (3) 

Here is the usual spacetime metric without the unimodular condition imposed. Let us 
consider a spatially compact spacetime M. On it we pick two nonintersecting spacelike 
slices of a spacetime history , which we may call E 2 and Si. These bound a region of 
spacetime 1Z. It follows that 



a- a = yf^g = Vof(K) (4) 
s 2 Jt, 2 Jk 

We may then call 

r(E) = fa (5) 

the elapsed spacetime four volume to the past of the surface E. 

The idea of OI6H3 is that the Hamiltonian constraints of general relativity are in this 
formulation replaced by evolution equations in this physically well defined time coordi- 
nate. To see explicitly what this means define the usual 3 + 1 slicing of spacetime, and set 
up a hamiltonian formulation of unimodular general relativity. In each surface E we can 
define a local time coordinate T = a*. This has conjugate momenta, tt, so that 

{f(x),n(y)} = fi 3 (y,x) (6) 

The hamiltonian constraints are then replaced by two sets of constraints. 

W = TVy/q -h = (7) 
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where h is the usual hamiltonian constraint, with vanishing cosmological constant, A, 
with density weight one. The second set of constraints are 

g a = d a 7r = o (8) 

This tells us that tt is a constant, and there are indeed equations of motion that confirm 
that it is a spacetime constant. Its value on any solution is the cosmological constant. 
From © we see that this variable cosmological constant is the momentum conjugate to 
elapsed four volume. 

In section 6 we discuss the implementation of (|7|8|) in an extension of the spin network 
representation. We show that the Wheeler-DeWitt equations (or quantum Hamiltonian 
constraint equations) can be interpreted to give a version of a many fingered time evolu- 
tion of the quantum state of geometry in the time r, in either the connection or the spin 
network representation. 

In the following we restrict ourselves to the case where spacetime has a topology £ x 
R where E is a compact three manifold without boundary. The case where there is a 
boundary or asymptotic conditions is interesting, but is reserved for a possible future 
paper. 



2 The Plebanski formulation of unimodular gravity 

Henneaux and Teitelboim[8] reformulated unimodular gravity so that the action depends 
on the full unconstrained metric and the gauge symmetry includes the full diffeomor- 
phism group of the manifold. They do this as follows. They introduce two auxiliary 
fields. The first is a three forrrfl a aj3l , whose field strength is b a ^ = da a ^s- The dual is a 
density^ 

b = jf^b^s = d a ~a a (9) 

where a a is the vector density field defined as a a = |e a/?7 ' 5 a / 3 7 5. The second is scalar field 
4> which serves as a lagrange multiplier. 

Their action translates into the Plebanski formulation as 

S HT = [ (B l A Fi - ^ ij B i A B j - <\)B % AB' + C matter ) + ^4>b (10) 
Jm 

where we impose the tracefree condition 

= (11) 

x We use greek indices for four dimensional manifold indices and latin indices for thee dimensional 
indices 

2 We use the notation where tildes refer to densities. 
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By varying cf) we find the unimodular condition emerging as an equation of motion 

b = Bi A B l (12) 

Varying a a we find 

a a = o (13) 

so that the field (p becomes a spacetime constant on solutions, so we can write 

<f>(x) = A (14) 
Varying B % we find the Einstein equations for any value of the constant A 

Fi = ABi + (15) 
plus matter terms. The last equation of motion is 

B i A B j - -5 ij B i AB k = (16) 
It makes sense to define the field 

% = ®H + \^<P (17) 
so that $ ' u = cj). If we use the equation of motion for B l we can rewrite the action as 

S HT2 = f ( F i A iTJ ($')-! + C matter ) + ^4>b (18) 

Jm o7rG 

Finally, we can give an interpretation of the field a as follows. Let us integrate ((12)) 
over a region of spacetime 1Z bounded by two spacelike surfaces Si and S 2 . Then we 
have 



a- a = Vol= y/^g (19) 

That is a pulled back into the surface is a time coordinate that measures the total four 
volume to the past of that surface. We can consider that time coordinate associated to a 
surface E to be 

t= J a. (20) 

3 Constrained hamiltonian dynamics of Plebanski unimod- 
ular gravity 

It is easy to construct the constrained hamiltonian framework for the Henneaux-Teiltelboim 
form of the theory (|T8|) . We do the usual 3 + 1 decomposition and define momenta for all 
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the fields. For this and the next section we neglect matter, it is trivial to reinsert it. We 
find the canonical momenta for the gauge field 

E? = c abc F bc (<$>% 1 (21) 

and also the primary constraints 

S = 7T - (f) = (22) 

as well as 

P$ = 11° = 7T = 7T C = (23) 

Here n a) n are the momenta conjugate to the a a and a , respectively. 
The definition ([22]) cannot be inverted to eliminate the velocities unless the following 
primary constraints are also imposed. 

V a = E\F\ h = (24) 

which generate diffeomorphisms and the modified Hamiltonian constraints 

U = e l > k E«E b F abk - <j)det{E«) = (25) 

We also find from the preservation of the vanishing of n° and tx c , respectively the two 
Gauss's law constraints 

G* = V a E ai = (26) 
G c = d c n = (27) 



where we have used (22). 



We can eliminate n and ti^ by using (|22|) with (|25|) to find 

W = e l]k EtE b 3 F abk - n det(E?) = (28) 
The standard calculation shows that there is also a Hamiltonian, which is 

ijk rpa rpb rp 

H = I (d a ~a«) I ^X: bk I (29) 



det(E? 

If one wants one can alternatively eliminate ((27|) and replace it by 

S c = d c \ — i — — = (30) 

V det ^ J 

Thus, the hamiltonian then vanishes on the constraint surface as ([30]) implies that there is 
for every solution a constant A for which 

e^ k EtE)F abk = k-det{Et). (31) 
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The Hamiltonian is then on the constraint surface equal to 

H = A J (d a a a ) « (32) 

The constraints rr a and G\ generate gauge transformations, respectfully 

~a a ^ a a + f a , a -> a + <9 c s c (33) 
where f a and 5 C are arbitrary vector densities. 

4 Physical observables of the classical theory 

We first show that there is in this formalism a definite method for calculating physical 
observables which is no more complicated than solving evolution equations as in uncon- 
strained dynamical theories. 

Let us eliminate the pair (a a , 7r a ) by solving the constraints n a = at the same time we 
choose a gauge generated by that constraint in which a a = 0. This gives us a formulation 
of the theory on the phase space Y defined by canonical pairs (A z a , E" 1 ), and (a°, n ). The 
constraints are W, (HHJ), V a (ED), Q 1 ^, Q c ((27} • Let us note that 

{r,7r (x)} = l (34) 

and write W as the undensitized form 

W = 7T - h = (35) 

where h is 

h = 1 2 (36) 

det(Ef) 

Let us now consider a physical observable on Y 

O = 0[A,E, a , 7v ] (37) 

Let us assume that it is locally gauge invariant and spatially diffeomorphism invariant so 
it satisfies 

{G\O} = (38) 
{V(v),O} = (39) 
for any vector field n" on E. We want to solve 

{W, O} = (40) 

{Ga,O}=0 (41) 
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We solve the latter, (|4~T]) first. This gives us 



d { dO 



dx a \da°(x) 



(42) 



This is solved by making O a function only of the integral r — f E a . Now we can consider 
(|40l> . This is one equation for every point x 6 E, 

= {W(x),0[t]} = ^ " {0[r\, h(x)} (43) 

However it is straightforward to show that if this is satisfied at any one point x G E it is 
satisfied for all points y e E. Consider that x and y are nearby so that there is a vector 
field v a such that h(y) = h(x) + {T>(v), h(x)} + where the ... indicate higher order terms 
in brackets with V{y) that come from exponentiation. Then we can write 

{0[r],h(y)} = {0[r],h(x)} + {0[T],{V(v),h(x)}} + ... 

= {0[r], h(x)} - {V(v), {h{x), 0[t)}} - {h(x), {0[t],V(v)}} 

= {0[r],h(x)} (44) 

where {D(v), {h(x),0[r]}} vanishes by d43j> and {0[t],V(v)} vanishes by p9). 

Thus we conclude that a complete set of conditions on 0[r] are the constraint equa- 
tions (|26|) , (j24]) together with the condition that it is a function of a° only through a depen- 



dence on the physical time, r together with the evolution equation 

dOW 



<)t {0[r],h(x)} (45) 



for any x. 



5 Construction of the path integral quantization 

The canonical theory we have arrived at has canonical pairs, (A l a , Ef), (a a , 7r a ) and (a , n ). 
The constraints are W, (|28]) , V a (|24)>, (y' (|26l> , ^ c (|27|) and 7r c . The ^ c count as one constraint 
because 

<9 [a £ c] = (46) 

Associated to these 11 constraints are 11 gauge fixing conditions. The first seven are 
standard from the Ashtekar formulation of general relativity, the four new ones gauge fix 
the a a and a° and will be discussed below. 

There is also a Hamiltonian which is ((29)) which is non-vanishing on the constraint 
surface. 
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We follow the standard construct the path integral from the the gauge fixed Hamilto- 
nian quantum dynamics. The partition function is 

Z = J dA i JE^daPda a dn a dn S(W)5(V a )8(g i )8(g c )S(7r c )8( gauge fixing)De* F p 

dt ^ ( E«A\ + 7r a° + ix c t - (d a ~a a ) f J J (47) 

The first thing we do is to integrate over 7r . This uses up the delta function in W, 
implying the substition 

( e ^E«E]F abk \ 
vtq ->■ 4 (48) 

in the rest of the integrand. This turns the delta functional in Q c into one in S c so the 
partition function is now 



Z = j dA i JE«dd Q da a d'n c 8(V a )8{g i )8{S c )5{'n c )8{ gauge fixing)Det FP 

/f ( ( E a E^ F \ \ 

* L p« + {-ikstT ) (a ° + 9 < &C) + " a v <49) 

We now introduce a vector density w a to exponentiate the constraint S a , and integrate 
over the n c . yielding 



Z = J dA*ji%da?da?dw c 6(V a )6(g i )6( gauge fixing) Det FP 

/f ( f €^^E a E^F \ \ 



We now may shift the definition of a a by 

a a ->• a a - w a (51) 

This eliminates any dependence on w a in the integrand so we trivially do the integral 
over dw a . At the same time we exponentiate the two remaining constraints, V a and Q % 
with la grange multipliers, respectively, N a , the lapse and A l . We then have 



Z = J dA i a dAidE?dN a dd°da c 8( gauge fixing) Det FP 

L J 

det(Ef) 



/f ( ( E a E^ F \ \ 

dt J I E«A\ + ( '-J** J (a + d c a c ) + N a V a + A& j (52) 
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We now can introduce the lapse N, which is an density of weight minus unity, by 
means of an insertion of unity, in the form 

1= fdN5(N- ~ a ° + d / aC ) (53) 
J V Det{E ai )J 

This gives us 

Z = J dA i a dAidEfdN a dNda°da c S (^N - g^t|^Q 5( gauge fixing)Det FP 

xexpi J dt J (E?Ai + Ne l ' k E1E b 3 F abk + N a V a + (54) 

We can change variables from E a \ N and N a back to the frame fields so we have 

Z = J dAidAide^daPdtfS (det(e) - a + d c h c ) 5( gauge fixing) J Det FP 

xexpi j dt J (e M A e u A F+) (55) 
5.1 Gauge fixing 

We now can introduce the four gauge fixing conditions 

fo = a°- t(e - d c a c ) = 0, f = a c = (56) 

where t is the time coordinate. e is a fixed density needed to make the density weights 
balence out. Then we have 

5° + d c a c = e . (57) 

The path integral now becomes 

Z = J dA^de^S (det(e) - e ) 8( gauge fixing)L>et FP 

xexpi J dt J (e M A e v A F+ ) (58) 

so we return to an action of the form of unimodular gravity. 

From here we can work back to the Plebanski formalism following the usual steps, 
giving us 

Z = / <MW^ (£, A B, - £0 ) ,(..),< gauge fixing) Drf < FP 

xexpijdtj (B i A Fi - ^ ij B i A B j ) (59) 

Thus, the path integral is the same as is usually taken to generate spin foam models, with 
the additional constraint that B l A B { is everywhere fixed to a background density e . 

10 



6 The unimodular loop representation 



Unimodular gravity in the Ashtekar variables was studied by 0, which we take as a 
starting point. Their work was done at a time when the roles of the volume operator and 
spin network basis were not fully appreciated, so it should be now possible to take their 
results some steps further. Here we just make a few remarks to motivate further work. 



6.1 The connection representation 

To quantize in the loop representation we may begin with an extended connection rep- 
resentation where the configuration space is a functional of connections A 1 on a three 
manifold S plus a scalar density field a . (As a a plays no role we will eliminate it and 
its conjugate momenta, which is a constraint.) For convenience we will rename aP, T and 
recall that r = f T measures total four volume since the beginning of the universe. 

Wavefunctionals are then functionals ty(A,T). The diffeomorphism and gauge con- 
straints affect only the A dependence in the usual way. in addition, there is the hamilto- 
nian constraint (|25l> which may be written, in singly densitized form 

d « z 

ih—,det{e)m(A, T) = hV(A, T) (60) 
dT 

where h is the singly densitized form of the hamiltonian constraint 

h = ^ k E«E)F ahk (61) 
There is the additional constraint Q c (|27)>, which may be written. 

ihd c -t^(A,f) = (62) 

This is solved by writing, 

*(A,f) = *(A,r(f)) (63) 
Alternatively we can separate variables, by writing 

*(A,f) = J d\(x)^(A,X)e l ^ fx{x) (64) 

which is a functional integral over values of a function X(x) on S. We then find that locally 

h(x)V(A, A) = \(x)det(e)V(A, A) (65) 
together with the condition that 

(W(i,A) = (66) 
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We see that *&(A, A) has support only on configurations where A is constant. 
We can integrate (|65|) over any region 1Z of the spatial manifold to find 



h n V(A, A) = \Vk*(A, A) (67) 

where 

h n = [ h (68) 
Jiz 

and 

Vk= [ Vq (69) 
If we define that constant value X(x) = A we find that the wavefunction evolves as 

V(A,t) = [ dA^(A,A)e lrA (70) 



This is now an ordinary integral over A. 

It is amusing to note that the Kodama state |12l| ^k(A, A) = e a / y cs{A) ^ g a so i u tion 
to (|65l> , with the state considered a function of variables A and A. With the standard point 
split regularizations, it solves the ordering 

e ijk k E) hahk ~ ^ a6 c4 c ) *(4 A) = (71) 

Whether this offers any improvement of the interpretational issues facing the Kodama 
state is unclear. 



6.2 The spin network representation 

We may now transform to the spin network representation, it is convenient to use the fact 



that in (|25|) , ir multiplies detE so we only need to have values of T on vertices of valence 
four or higher, on which the volume operator doesn't vanish. Let TZi be a decomposition 
of S into regions, each enclosing a single vertex Vi of a spin network state of valence four 
or higher. Then we can associate to each vertex v it the integral 

f = n (72) 

so that Yl,i T i = T which measures the past four-volume of the slice. Note that the par- 
tition of r so defined on vertex labels is arbitrary, as the regions are arbitrary parts of a 
regularization procedure. 

We then augment the usual definition of spin network states with these labels 7j on 
each vertex of non-vanishing volume. The quantum states are then functionals 

(73) 
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But the solutions to (|27) impose that all the 7$ = r so we have quantum states of the form 

(74) 

that satisfy for every region 1Z, 



ih^v 1l *(r,T) = h R *(r,T) (75) 



7 Conclusions 



In this paper we have extended to the Ashtekar and Plebanski formulations of gravity 
the results of llTOl , in which it is shown that the quantum effective action of unimodular 
gravity retains the unimodular property, ©. The main results of this paper are the formal 
expression for the path integral in unimodular form of the Plebanski action, eq, (|59l> and 
the unimodular form of the Wheeler-deWitt equation in spin network representation, eq. 



§75) . We note that the when the region 1Z in the latter is chosen to be the whole spatial 
manifold, we have a kind of Schroedinger equation. But we should caution that this 
equation has to be satisfied for every region, 1Z, of the manifold. In particular, note that 
there are generically regions in spin network states in the kernel of the volume operator 
where the regularized forms of hn so far studied act non-trivially. These include nodes of 
valence three; on these the hamiltonian constraint will take its usual form. What is new is 
the action on regions, such as those containing generic nodes of valence four or higher, on 
which the action of both Vr and hn is non-trivial. The development of solutions in these 
regions will give something like a many fingered time, where the "extent" of each finger 
is determined by the action of the volume operator. 

Hence, while unimodular gravity does offer a physical time coordinate, evolution in 
terms of which is defined by a Schroedinger equation, this does not get us out of the 
challenge of solving an infinite number of dynamical equations, one for each point or 
region of the spatial manifold. 

Given that there has been recently a great deal of progress in the construction and 
analysis of spin foam models, it would be very interesting to see if a concrete definition 
can be given in spin foam terms of the path integral (|59]) . We note that some progress 
in this direction has been made in |TT|| which also studies the unimodular form of loop 
quantum cosmology. 
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